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ABSTRACT 

We attribute the observed rotational velocity profiles of spiral galaxy disks to 
a combination of Newtonian and exponential gravitational potentials. We offer 
a novel, yet perfectly plausible, Lagrangian from which the exponential potential 
is derived. The exponential potential is defined by two universal constants (cor- 
responding to a coupling constant and a graviton mass) that we determine from 
our sample of 12 THINGS disks. All velocity profile fits are excellent. 

Subject headings: dark matter — elementary particles — galaxies: kinematics 
and dynamics — galaxies: spiral — gravitation 

1. Introduction 

Kepler's third law states that the squares of the periods of any two planets are in the 
same proportion as the cubes of their mean distances from the Sun; this is explained by the 
Newtonian inverse-square law of gravitation. This relationship does not apply to the stars 
in a spiral galaxy disk because their periods are generally more or less proportional to their 
distances from the galactic center. That is, the disk circular velocity Vc profiles of spiral 
galaxies are observed to be approximately flat, Vc ~ constant. Newtonian gravitation does 
not explain this, so astrophysicists have been disposed (1) to hypothesize the existence of 
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invisible matter in a halo centered on each galaxy, and that this dark matter is distributed 
in a seeming conspiracy to result in the observed motions while retaining the inverse- square 
law of gravitation, or (2) to modify Newtonian theory. 



The Lagrangian for Newtonian gravitation in a centrally symmetric field is 



C 



where G is the gravitational constant, ip = i^i^) is the gravitational potential, and p = p{r) 
is the mass density. The solution to the variation 5£ = is the Poisson equation. 



V^^ = AnGp, 



(2) 



so the solution to Eq. [2]for a point mass M at r = is the Newtonian p otential V ' y = —G M/ r, 
and the consequent centripetal force obeys the inverse-square law. Milgrom ( 1983al Jbl) ad- 
dressed the disk velocity discrepancy for gala xies with his Modified Newtoii ian Dynamics 
theory, MOND, for which Eq. [T] is modified (IBekenstein and MilgromI 1 19841 ) with the re- 
placement 

(V^)2^a2j-[(V^/ao)2], (3) 

where the acceleration is a universal constant and ~ x for x ^ 1 and ~ 1 

for X ^ 1. On setting Oq ~ 10~* cm s~^, the functional form of J-" o bliges spiral galaxy disk 



velocity pr ofiles to be generally flat, i n agr eement with observations. iBrownstein and Moffat 



(120061 ) and iMannheim and O'BrienI (120111 ) also devised non- Newtonian explanations for the 
fiat velocity profiles of spiral galaxies. 

Eckhardt. Pestana fc Fischbach 1 J2010h [EPF] considered non-Newtonian gravitation at 



the scale of galaxy superclusters. For the Lagrangian, 



(4) 



the solution to 6C = is 



(V2-/x2)^ = 47rGp, (5) 

and the solution to Eq. [5] for a point mass M at r = has the form of a Yukawa potential 
(see Appendix), 

ipY = —GMexp{—pr)/r. (6) 

If /i > 0, this results in a Milne universe that is in full accord w ith cosmological ex pansion ob- 
servations of type la supernovae but, unlike the explication of iRiess et al.l (120071 ). it requires 
neither dark matter [Qrn = 0) nor dark energy (^a = 0). The choice ~ 5 Mpc then 
explains the scales of galaxy superclusters and of the fundamental spectrum of the cosmic 
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background radiation; and it explains why neighboring superclusters tend to be ahgned in 
spongiform "surfaces" with vast empty regions between them. 

The graviton mass corresponding to fi^^ = 5 Mpc is m = hfi/c = 1.3 x 10^'^° eV/c^, 
so EPF conjectured that the mass of the (cosmological) graviton is ~ 10^^*^ eV/c^, and 
that there could be heavier gravitons as well, but none that is lighter than rric. This led us 
to conjecture that a heavier (galaxy) rUg graviton is responsible for t he flat galaxy vel ocity 



curves. An rUg Yukawa potential cannot explain the flat curves, but lEckhardtl (119931 ) had 
suggested that an rUg exponential potential combined with a Newtonian potential could 
provide an explanation (see his Fig. 2). However, his theoretical argument for the existence 
of an exponential potential was arduous and unconvincing, so we now offer a more substantial 
argument based on the Lagrangian, 



(7) 



The solution to 5£ = is 

/i-^VV - 2VV + = Ai~^(V^ - /i^) V = -47rGp, (8) 

and the solution to Eq. |8]for a point mass M at r = has the form of a exponential potential 
(see Appendix), 

V'e = — T^M/i exp(— yur), (9) 
where, for the irig graviton, fi ^ fig = nigc/h and 7 — > 7g is a dimensionless constant. 



2. Modeling Galaxy Disks 

At galaxy scales, Yukawa potential due to appears to be Newtonian, so we examine 
how well the combination of Newtonian and rrig exponential potentials can reconcile obser- 
vations of luminosity distributions of spiral galaxies with those of the rotational velocity 
distributions of the gases and stars in their disks. 

Because the centripetal acceleration of a disk object rotating in a circular orbit at 
velocity Vc and distance r from a galaxy centroid equals the specific force on the object, 

v'^ = ri^'j, + ij'j,) . (10) 

We model the Newto n ian a nd nig potentials of the disk using the double ring technique of 



Eckhardt &: Pest ana ( 2OO2I ) [EP], using their Eq. 3 and performing the calculations with 



Mathematica, Version 8.0. Adopting the EP [c„, R, a] nomenclature, we set c_i = 1 to 
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evaluate ipN, whereas to evaluate ipE we use the series expansion of exp{—figR) to calculate 
c„ for n = 0, 1, • • ■ , 10; and we then scale up the final result by the factor jg. The density 
distribution (j{R) is calculated using all photometrically determined galaxy matter in its 
central component and disk, plus interstellar gases in the disk, which are determined from 
the 21-cm radiation of neutral hydrogen - along with a 25 percent by mass admixture of 
helium. Our sample of spiral galaxies is a subset of 12 nearby disks selected from the 19 



disks in the HI Nearby Galaxy Survey [THINGS] that were analyzed by lde Blok et al\ (120081 ) 
[dBETAL]. We excluded seven galaxies a priori either because they are clearly dominated 
by noncircular motions (NGC 3031, NGC 3627, NGC 4826 and NGC 4736) or we cannot 
eliminate that possibility (NGC 2366), or because of the the difficulty and uncertainty of 
their photometric analyses (NGC 3521 and NGC 3198). 

The galaxies in our sample, which were all observed with the same instrument and 
resolution, extend over a luminosity range of more than three orders of magnitude. They 
comprise the best currently available observational data for resolving galaxy rotation profiles 
that can be used to evaluate diverse gravitational models. dBETAL performed a meticulous 
decomposition of the stellar surface density profiles that is of utmost importance for testing 
non-Newtonian dynamical theories; indeed, they even weighted each point of the observed 
photometric profile by the M/L ratio expected from the observed color gradient as a function 
of radius. Table [1] lists our sample in order of decreasing rotation velocity. Bear in mind 
that these astronomical data, although highly refined, are not perfect: there always remain 
lingering uncertainties which, in order of decreasing significance, arise from errors in distance, 
photometric imprecisions, and sub-estimates of radio observations. 

To reconcile the velocities, Vc, determined by Eq. [10] with velocities, Vobs, observed using 
the 21-cm lines, we use the observed surface brightness data and color gradients without any 
adjustments, and the HI densities of dBETAL by applying the double ring technique. We 
use this technique for the principal disk and the inner disks (bulges) detected by dBETAL. 
Then we iteratively adjust the mass-luminosity ratios, M/L for the Spitzer IRAC 3.6 /im 
band, in each of the 12 dBETAL galaxies' disks, and central components where appropriate, 
along with the two universal constants, 7^ and fig, until the Vc and Vobs profiles are in close 
agreement for all galaxies. Note that the distances are fixed parameters (the same as those 
used by dBETAL for their dark halo models) which have uncertainties of > 10%. Allowing 
these distances to vary within those uncertainties would result in modest improvements to 
the already excellent fits. We determine that the numerical values of the two fundamental 
constants are 7^ = 12.5 and fij^ = 20 kpc (m^ = 3.2 x 10^^*^ eV/c^). Although we have 
no formal error estimates for these constants, they are narrowly constrained. Table [2] lists 
our dynamically determined values for M/L. There are only one or two free parameters per 
galaxy, depending on whether (Figured]) or not (Figure ]2]) the galaxy has an additional inter- 
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nal photometric disk. By contrast, the dark matter hypothesis requires three free parameters 
per galaxy, or four parameters for the galaxies with inner disks. 

The quality of the adjustment can be assessed by inspecting the curves of Figures [T] 
and [21 and by the plausibility of the dynamically determined ratios M/L. The velocity 
components in Figure [1] and |2] combine as root-sum-squares. Each of the 12 panels of the 
figures shows that the rrig potential component makes up for the discrepancy between the 
Newtonian component and the observed velocities. The velocity Vc profiles agree very well 
with the observed velocity Vobs profiles, for high as well as low Vobs, and especially so in 
view of observational complications and the simplicity of our model. Also, our fits indicate 
that there is a smooth transition between the small and large disks, with the mg potential 
becoming progressively more important as the radial dimensions of the galaxies increase. 



or 



For compariso ns with Figures [T] and 121 dark matter halo fits are shown by dBETAL 



Oh et all feOOSh and, exc e pt fo r NGC 925 and IC 2574, the MOND fits are shown by 



Gentile. Famaey &: de Blok I (120111 ). Note that each dark matter halo fit requires a pair 
of parameters that is idiosyncratic to the particular galaxy being modeled, 2 x 12 = 24 
parameters in all, and that the MOND fits require recurrent tweakings of both the "universal 
constant" Qq and the galaxy distances (iBottema et al.ll2002[ ). By contrast, all of our fits are 
for the same two constants, jg and fig, and the galaxy distances are not modified. 

All the dynamically determined M/L ratios in Table |2] are unambiguous and, with one 
exception (DDO 154), they are consistent with stellar population synthesis models. The same 
exception was noted by dBETAL in their dark halo modeling; we refer the reader to their 
paper for a discussion of possible explanations. Also, the inner M/L ratios for NGC 7331 and 
NGC 2903 appear odd. However, possible explanations could be the presence of a strong dust 
ring in the inner disk of NGC 7331, and the presence of the ba r and a very dense molecular 



disk in the inner disk of NGC 2903. In fact, both dBETAL and lGentile. Famaey fc de Blok 
( 1201 ll ) also get unrealistic M/L ratios for these two galaxies. 



3. Conclusions and Outlook 

The essential result of this work is that the mg exponential potential can account for 
the magnitude of the discrepancy and reproduces the general shapes of the rotation curves 
of these 12 THINGS spiral galaxies. Invoking dark matter halos to reproduce the observed 
shapes requires 24 additional free parameters which have no underlying physical significance; 
they are chosen just to get satisfactory fits. And, although MOND purportedly involves only 
one universal constant, ao, it is actually a parameter that does not remain precisely the same 
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Table 1. Model Inputs. Luminosities are in the 3.6 nm band. 



Object 


Distance 
[Mpc] 


Disk Luminosity 
[10' L©] 


Inner Disk Luminosity 
[10' Le] 


HI Mass 
[10' Mq] 


NGC 2841 


14.10 


173.15 


21.54 


98.91 


NGC 7331 


14.72 


398.14 


1.66 


75.81 


NGC 6946 


5.90 


92.61 


0.02 


37.26 


iNGG 2903 


8.90 


30.34 


1.87 


4.63 


NGC 5055 


10.10 


169.25 


19.77 


98.9 


NGC 3621 


6.64 


39.76 




67.44 


NGC 2403 


3.47 


16.46 


0.83 


30.42 


NGC 925 


9.20 


17.12 




28.37 


NGC 7793 


3.91 


9.41 




7.99 


NGC 2976 


3.56 


4.60 




1.13 


IC 2574 


4.00 


2.24 




13.32 


DDO 154 


4.30 


0.82 




3.15 



Table 2. Model Outputs (in solar units). 



Object 


(M/L) 


(M/L)in„er 


NGC 2841 


0.10 


3.00 


NGC 7331 


0.02 


23.00 


NGC 6946 


0.34 


0.01 


NGC 2903 


0.15 


11.00 


NGC 5055 


0.06 


1.00 


NGC 3621 


0.25 




NGC 2403 


0.50 


2.50 


NGC 925 


0.40 




NGC 7793 


0.70 




NGC 2976 


0.70 




IC 2574 


0.70 




DDO 154 


9.50 
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Fig. 1. — Velocity fits for galaxies witli two photometric disks. Tlie tliick solid black line is 
V. The thin solid lines are the v contributions of the Newtonian potential and the broken 
lines are those of of the rrig potential. 
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Fig. 2. — Velocity fits for galaxies with just one photometric disk. 
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for all fits, and the galaxy distances and the exact functional forms of T (see Eq. |3]) are also 
sometimes modified. By contrast, the vrig exponential potential is determined by only two 
fixed and physically significant parameters: /x^, which is proportional to the exchange boson 
mass, and 7(,, which is proportional to the square of the coupling constant. The inequality 
2 < 24 provides a compelling argument for the superiority of the vrig exponential potential 
over ad hoc dark matter halos for modeling spiral galaxy rotation profiles. A similar, although 
not so lopsided, inequality apphes on contrasting the potential fits with those of MOND 

The Yukawa potential is derived from the Eq. HlLagrangian, and the rrig exponential 
potential is derived from the Eq. [7]Lagrangian. In the solar system, the mc potential appears 
to be Newtonian whereas the nig potential results in a Rindler type acceleration toward the 
Sun, 'fgfi'^GMQ = 4 X 10~ ^^ m s~^, w hich is far below the celestial mechanics detection 



threshold, 1 x 10"^^ m s'^ fllorioll201l[ ). 



A. Solutions to Eqs. [5] and [8] 

Modify the Fourier integral pair 

POO 

$(A;) = / (f){r) sinkr dr, 
Jo 



oo 



2 

(r) = — / sinkrdk 



TT 

to 



k Jq r 

jo{kr) k dk, 



(r) _ 2 r m ^ , , 



r n Jq k 
where the spherical Bessel function jo{kr) is simply 

, sin kr 

Next make the replacements ^{k)/k ^(fc) and (f){r)/r — )■ iplr). This gives a Fourier-Bessel 
integral pair: 

POO 

^(fc) = T[^(r)] = / tlj{r)jo{kr)r'^dr, (Al) 



2 /*oo roo 

il){r) = T-^[m{k)] = - j ^!{k)jQ{kr)k'^dk = — j "^{k) expikr kdk. (A2) 
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Then 
and 

T[V'^ilj{r)] = -P^(A;). 

Suppose that p = p{r) = po for r < tq, where kr^ <^ 1 (so jo{kr) ^ 1), and that p(r) = for 
r > ro; then 

pro 

T[AnGp] ^ Gpo / Ajrr^ dr = GM. 



Jo 

This becomes an equahty if M is a point mass, and 

T[(V^ - p^)iIj - AnGp] = -{k^ + /i^)^(A;) - GM = 0, 

or 

(A3) 

Contour integration in the upper half plane about the pole k = ip then gives the Yukawa 
potential of Eq. ^ ip = ipy- Following the same approach to solve Eq. [HI the only modification 
is that Eq. I A3 1 is replaced by 

p'^GM , 

m-'j^- (A4) 

Rather than perform a contour integration about a double pole, it is simpler just to note 
that 

-p^ ^ = Mk) 

so the corresponding potential is modified to 

7 2 ^V'y A* di^Y GM 

But G here is a generic constant anyway, with the same dimensions as the gravitational 
constant G used in Eq. Ej so we simply replace G/2 by 'yG, where 7 is a dimensionless 
constant that takes into account the 1/2 factor and the differences between the coupling 
constants for the different potentials. The solution to Eq. [8]is then the exponential potential 
given by Eq. [9l The p is not subsumed into 7 because p has the dimension of inverse distance, 
as does 1/r in the Yukawa potential. 
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